Bayesian Inversion of Time Domain Electromagnetic (TDEM) Synthetic Data using Langevin Monte Carlo (LMC)
نویسندگان
چکیده
The inversion of time domain electromagnetic (TDEM) data is an ill-posed problem. This makes the standard procedure, which to linearize related objective function. Then take a deterministic approach determine solution that can minimize said function, has potential be trapped in local minimum. In this study, we solved problem TDEM using Bayesian framework by generating sample from posterior distribution. distribution contains information uncertainty and results forward modeling formulation, prior about subsurface parameter model. conducting sampling process, use Langevin Monte Carlo (LMC) algorithm, one many gradient-based Markov chain (MCMC) algorithms. was performed on synthetic generated through several test models. We used model with varying thickness resistivity values for inversion. also added probability smoothness constraint between adjacent layers allowing sharp smooth transitions.
منابع مشابه
Transdimensional Monte Carlo Inversion of AEM Data
The majority of existing methods used for the inversion of airborne electromagnetic (AEM) data use what are generally called gradient-based optimization techniques. They typically minimize an objective function comprised of data misfit (e.g. least squares) and model regularization (e.g. roughness) terms. Since the problem is non-linear, an iterative search involving the matrix solution of equat...
متن کاملMonte Carlo Noise Reduction Using Bayesian Method in Wavelet domain
A novel post-processing approach for removing Monte Carlo noises in synthetic images is presented in this paper. This paper first presents our findings on the statistical characteristics of the Monte Carlo noise, and then proposes a Bayesian method to remove this noise. The aim of this approach is to efficiently produce high quality synthetic images using Monte Carlo based rendering at low samp...
متن کاملFinite-Time Analysis of Projected Langevin Monte Carlo
We analyze the projected Langevin Monte Carlo (LMC) algorithm, a close cousin of projected Stochastic Gradient Descent (SGD). We show that LMC allows to sample in polynomial time from a posterior distribution restricted to a convex body and with concave log-likelihood. This gives the first Markov chain to sample from a log-concave distribution with a first-order oracle, as the existing chains w...
متن کاملRiemann Manifold Langevin and Hamiltonian Monte Carlo
This paper proposes Metropolis adjusted Langevin and Hamiltonian Monte Carlo sampling methods defined on the Riemann manifold to resolve the shortcomings of existing Monte Carlo algorithms when sampling from target densities that may be high dimensional and exhibit strong correlations. The methods provide fully automated adaptation mechanisms that circumvent the costly pilot runs required to tu...
متن کاملStochastic Quasi-Newton Langevin Monte Carlo
Recently, Stochastic Gradient Markov Chain Monte Carlo (SG-MCMC) methods have been proposed for scaling up Monte Carlo computations to large data problems. Whilst these approaches have proven useful in many applications, vanilla SG-MCMC might suffer from poor mixing rates when random variables exhibit strong couplings under the target densities or big scale differences. In this study, we propos...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Journal of physics
سال: 2022
ISSN: ['0022-3700', '1747-3721', '0368-3508', '1747-3713']
DOI: https://doi.org/10.1088/1742-6596/2377/1/012023